Brownian transport in corrugated channels with inertia 
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The transport of suspended Brownian particles dc-driven along corrugated narrow channels is 
numerically investigated in the regime of finite damping. We show that inertial corrections cannot 
be neglected as long as the width of the channel bottlenecks is smaller than an appropriate particle 
diffusion length, which depends on the the channel corrugation and the drive intensity. Being such 
a diffusion length inversely proportional to the damping constant, transport through sufficiently 
narrow obstructions turns out to be always sensitive to the viscosity of the suspension fluid. The 
inertia corrections to the transport quantifiers, mobility and diffusivity, markedly differ for smoothly 
and sharply corrugated channels. 

PACS numbers: 05.40.-a,05.60.Cd,51.20.+d 



I. INTRODUCTION 

Brownian transport in narrow corrugated channels is a 
topic of potential applications to both natural [IH2] and 
artificial devices 0]. Depending on the amplitude and 
geometry of the wall modulation, corrugated channels 
fall within two distinct categories [Fig. [l] : (i) smoothly- 
corrugated channels (e.g., as shown in Fig. 1(a)). Also 
called entropic channels [5], these quasi-one-dimensional 
(ID) channels were first introduced in Ref. [6 and further 
investigated in Refs. [7rH3]. as an instance of two (2D) 
or three dimensional (3D) systems describable in terms 
of an effective ID kinetic equation. These are typically 
modeled as periodic channels with axial symmetry and 
unit cells delimited by bottlenecks which are assumed 
to be narrow with respect to the cell dimensions, i.e., 
the channel cross section and the period; (ii) compart- 
mentalized (or septate) channels [T4lll9| . These septate 
channels are sharply corrugated channels formed by iden- 
tical compartments separated by thin dividing walls and 
connected by narrow openings (pores) centered around 
their axis. At variance with the case of smoothly cor- 
rugated channels, diffusion in compartmentalized chan- 
nels cannot be reduced to an effective ID kinetic pro- 
cess directed along the axis. Accordingly, driven trans- 
port in such strongly constrained geometries exhibits 
distinct features, which cannot be reconciled with the 
known properties of Brownian motion in quasi- ID sys- 
tems |3 1201 12] • 

Corrugated channels are often used to model transport 
of dilute mixtures of small particles (like biomolecules, 
colloids or magnetic vortices) in confined geometries [3]. 
Each particle is subjected to thermal fluctuations with 
temperature T and large viscous damping constant 7, 
and a homogeneous constant force directed locally paral- 
lel to the channel axis. Such a dc drive is applied from the 
outside by coupling the particle to an external field (for 
instance, by attaching a dielectric or magnetic dipole, or 
a magnetic flux to the particle), without inducing drag 
effects on the suspension fluid. Interparticle and hydro- 
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FIG. 1: (Color online) Sketch of a smoothly corrugated (a) 
and a compartmentalized (b) 2D channel directed along the 
a;-axis. In both cases the channel unit cell is xl long and yi 
wide; the radius of the connecting bottlenecks or pores is A. 



dynamic interactions are thus ignored for simplicity (for 
a discussion on this assumption see Ref. [5]). 

In this paper we investigate the relevance of the iner- 
tia effects due to the viscosity of the suspended particle. 
As is often the case with biological and most artificial 
suspensions [4], the Brownian particle dynamics in the 
bulk can be regarded as overdamped. This corresponds 
to (i) formally setting the mass of the particle to zero, 
m — 0, or, equivalently, to make the friction strength 
7 tend to infinity, and (ii) assuming F smaller than the 
thermal force Fq = "f^kT/m (Smoluchowski approxima- 
tion) |22j . The current literature on corrugated channels 
invariably assumes such an overdamped limit. But how 
large is an infinite 7 (or how small a zero m)? The an- 
swer, of course, depends on the geometry of the channel. 

Our main conclusion is that the overdamped dynam- 
ics assumption for Brownian diffusion through pores of 
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FIG. 2: (Color online) Rescaled mobility, 7/1, in a smoothly corrugated channel with r = xl/vl = 1, A/2/1, = 0.1, and (a) 
versus 7/7F for different J* 7 ; (b) versus F/ Ft for different 7. The relevant scaling parameters are Ft = kT/A, jt = VmkT/ A, 
Eq. |8|), and 7f = ^/ mF/A, Eq. ([9]). The dashed lines represent, respectively, the fitting power laws (7/7f) q in (a) and 
(F / Ft) - "'' 2 in (b), both with a = 1.4 (see Sec. Ill I. In (c) 7/i is plotted versus 7/7F (main panel) and F/Fq (inset) for 
0.8, F/Ft = 125, and different cross-section ratios, A/yt- The corresponding fitting exponents a are also reported 
In the inset, F is expressed in units of Fa instead of Ft for graphical reasons. The dependence of 7/1 on the 



in the legend. 

geometry of the channel unit cell for low damping and small drives is illustrated in ( d), w here 7/1 is plotted versus r*y/~/T- The 
predicted linear law with slope n/4 [55] is represented by a dotted line [see also Eq. (12 1 and text following]. 



width A subjected to a homogeneous drive F, applies 
only for the conditions [22] : 

7 > VmkT /A, and 7 > ^mF/A, 

irrespective of the degree of corrugation. This means 
that inertial correction cannot be neglected as long as 
the Brownian diffusion is spatially correlated on a length 
(It = VrnkT/j at small dc drive, or lp = mF/7 2 at large 
dc drive) of the order of, or larger than, the pore width 
A. Therefore, for sufficiently narrow pores or sufficiently 
large drives, inertia always comes into play by enhancing 
the blocking action of the channel bottlenecks. 

This paper is organized as follows. In Sec. [TT] we in- 
troduce the Langevin equation formalism employed in 



our simulation code. Simulation data for the particle 
mobility and diffusivity are analyzed in Sec. |III| as func- 
tions of the drive, the channel geometry, and the damping 
constant in sinusoidally-corrugated channels. We report 
significant deviations from the best known overdamped 



regime. In Sec. IV we consider the case of septate chan- 
nels for which dependable fitting formulas could be an- 
alytically obtained. Inertial effects in these two limiting 
corrugation regimes are compared in Sec. |Vj Finally, in 



Sec. VI we add some concluding remarks. 
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II. MODEL 



diffusivity, 



Let us consider a point-like Brownian particle of mass 
m diffusing in a 2D suspension fluid contained in a pe- 
riodic channel with unit cell xl x yL, as illustrated in 
Fig. [T] The particle is subjected to a homogeneous force 
F. The damped dynamics of the particle is modeled by 
the 2D Langevin equation, 



l dt 2 



dr 



J—+F +v^fcT|*(t), 



dt 



(1) 



where r — (x,y). The random forces = (£,x{t),£, y (t)) 
are zero-mean, white Gaussian noises with autocorrela- 
tion functions (^i(t)^j(t')} = 2SijS(t — t'), with i,j = x,y. 
Here, 7 plays the role of an effective viscous damping con- 
stant incorporating all additional effects that are not ex- 
plicitly accounted for in Eq. ([!]) , like hydrodynamic drag, 
particle-wall interactions, etc. We numerically integrated 
Eq. ([I]) by a Milstein algorithm . The stochastic aver- 
ages reported in the forthcoming sections were obtained 
as ensemble averages over 10 6 trajectories with random 
initial conditions; transient effects were estimated and 
subtracted. 

As anticipated in Sec. |TJ we considered two categories 
of periodic channels, smoothly-corrugated and septate 
channels. The symmetric walls of smoothly-corrugated 
channels have been modeled by the sinusoidal functions 
±w(x), where 



w\x) 



{y L + A) - (y L - A) cos 



/2ttx\ 



(2) 



as shown in Fig. [fja)]. On the other, the compartments 
of the septate channels are rectangular and the dividing 
walls have zero width, as shown in Fig.[]Jb). 

Two quantifiers have been used to best represent the 
different transport properties of these two channel ge- 
ometries in the overdamped limit, 7 — > 00: 

(i) mobility. The response of a Brownian particle in 
a channel subjected to a dc drive, F, oriented along the 
x-axis is expressed by the mobility, 



H(F) = (v(F))/F, 



(3) 



where (v) = (x{F)) = lim f ^ oc [(a;(<)) - x(0)]/t. In en- 
tropic channels, /i(F) increases from a relatively small 
value no, for F — 0, up to the free-particle limit, 7/^00 = 
1. for F — > 00 [IT]. We recall that in a smooth chan- 
nel a free particle drifts with speed = F/7, that is, 
with 7/1 = f. On the contrary, in compartmentalized 
channels jfi(F) decreases monotonically with increas- 
ing F towards a geometry-dependent asymptotic value, 
7/ioo, equal to the ratio of the pore to the channel cross- 
section [TJ1, that is 



7^00 = A/y L , 



(4) 



(ii) diffusivity. As a Brownian particle is driven across 
a periodic array of bottlenecks or compartment pores, its 



D(F) = lim 

t—>CO 



(x 2 (t)) 



(x(t)) 2 ]/2t, 



(5) 



picks up a distinct F-dependence. In entropic channels 
with smooth bottlenecks, for F — » 00 the function D(F) 
approaches the free diffusion limit, -D(oo) = Dq, after go- 
ing through an excess diffusion peak centered around an 
intermediate (temperature dependent value of the 
drive. The bulk or free diffusivity, D , is proportional 
to the temperature, Do = kT/mj. Such a peak signals 
the depinning of the particle from the entropic barrier ar- 
ray |24j . In compartmentalized channels, instead, D(F) 
exhibits a distinct quadratic dependence on F [T5J [T7] , 
reminiscent of Taylor's diffusion in hydrodynamics [25] . 
that is, for A <C Ul, 



D(F) 



I /FA 
F» 2 I ~kf 



(6) 



This observation suggests that the particle never frees 
itself from the geometric constriction of the compartment 
pores, no matter how strong F. 

These two quantifiers can also be used to assess the 
magnitude of the inertia effects. We remind here that, 
in the absence of external drives and for any value of the 
damping constant, Einstein's relation |21j 



7^0 = D(0)/D o , 



(7) 



establishes the dependence of the transport parameters 
on the temperature and the channel compartment geom- 
etry under equilibrium conditions [14] . 

In preparation for the quantitative analysis of our nu- 
merical data, we remark that Eq. can be conve- 
niently rewritten in terms of the rcscalcd units t — > jt/m 
and x/It, with It — %/ mfcT/7. A straightforward di- 
mensional argument shows that, for any given chan- 
nel unit cell xl x y^, both the particle rescaled mo- 
bility, 7/!, Eq. pj, and its rescaled diffusivity, D/D , 
Eq. ([5]), are functions of the rescaled drive, F/F , with 
Fo = ^jyJkT/m, and three cell parameters, typically, the 
pore width, A/ It, the pore-to-channel cross-section ratio, 
A/yL, and the compartment aspect-ratio, r = x^/yi, (see 
Table [T]). Note that a simultaneous rescaling of all lengths 
by a factor k would correspond to a noise intensity rescal- 
ing, T —> T/k 2 . Throughout our simulations we assumed 
narrow channels with small bottlenecks, meaning that 
xl > Ul and A < y L . 



III. CORRUGATED CHANNELS 

As anticipated in a preliminary report [26] , inertial ef- 
fects in corrugated channels become apparent both for 
small 7 and for large F. By inspecting Fig. [2] we imme- 
diately realize that inertia tends to suppress the particle 
mobility through the channel bottlenecks. Indeed, in the 
underdamped limit, 7—5-0, the rescaled mobility drops 
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TABLE I: Summary of characteristic scaling parameters and their meaning 



ro = 7a/ — 
1 y m 


Thermal force: viscous force experienced by a Brownian particle 
with thermal velocity Vth = \JkTjm. 


n kT F 

L>a = — , Woo = - 


Free diffusivity and velocity in bulk. 


; VmfcT 
IT = 

7 


Thermal length: distance covered by a Brownian particle diffusing 
with thermal velocity Vth in the relaxation time, m/7. 


; mF 
t F — — T 


Ballistic length: distance covered by a driven Brownian particle 
drifting with velocity «oo in the relaxation time, m/7. 


VmkT 

7T = A 


Damping cut-off at the pore (zero drive): It — A. 


/mF 
7 p = V - A" 


Damping cut-off at the pore (strong drive): If = A. 


= f , Ft = % 


Scaling parameters introduced in Figs. 2|5 

obtained by replacing 7 with jt, respectively, in Do and Fq. 



to zero, no matter what F [Fig. [2ja)]. In particular, 
when expressing 7 in units of jp — yj mF/ A, see Eq. ^ 
below, the mobility curves at large drives tend to collapse 
on a universal curve well fitted by the power law (7 /7f) q 
with a = 1.4. Correspondingly, in Fig. (2^b) the mobility 
decays like F~ a / 2 for small F >• Ft- 

The power law, 7/1 oc (7/7f) q , introduced here is only 
a convenient fit of the rescaled mobility function, even if 
it holds for two or more decades of "f/jF- [Note that the 
power law 7/i oc F~ a / 2 works throughout the entire F 
range explored in Fig. [2jb) .] The analytical form of that 
function remains to be determined. The data reported 
in Fig. [2jc) clearly suggests that the fitting exponent a 
slightly depends on A, with a — > 2 in the limit A — > 0. 

The dependence of the rescaled mobility on the sys- 
tem parameters in the underdamped limit is further il- 
lustrated in Fig. J2jd) , where at low 7 and for vanishingly 
small drives, the mobility grows proportional to the as- 
pect ratio r = xl/i/l. of the channel unit cell and the 
pore cross section A. 

Deviations from the expected overdamped behavior are 
the more prominent in the diffusivity data. As shown in 
Fig. [3j at large 7 the curves D(F) approach the horizon- 
tal asymptote D(F) = D , as expected [TT]. However, 
beyond a certain value of F, seemingly proportional to 
7 2 (inset), these curves abruptly part from their horizon- 
tal asymptote with a sort of cusp. In the underdamped 
limit, the F dependence of the diffusivity bears no resem- 
blance with the typical overdamped behavior. At low 7, 
all D(F) data sets collapse on a unique curve [Fig. |3j in- 
set], which tends to a value smaller than Dq for F — >• 0, 



and diverges for F — > 00, like F^ with j3 ~ 1. Such power 
law holds for large 7, as well, though for sufficiently large 
F, only. Indeed, for exceedingly large F, all D(F) curves 
seem to eventually approach a unique asymptote, irre- 
spective of 7. 

By comparing the plots of Figs. [2][3] we conjecture 
that corrections due to inertia become significant in two 
regimes, namely: 

(i) at low drives for 

7 < 7 T = VmfcT/A . (8) 

This characteristic damping was used to rescale the mo- 
bility data in Fig. [l] [see also Fig. [2jb), inset]; moreover, 
in Fig. [3j for 7 < 7t the diffusivity becomes a monotonic 
function of F with no plateau around Dq. The physi- 
cal meaning of 7t is simple. For 7 < 7t the thermal 
length It = V mkT/j grows larger than the width of the 
pores, A, so that the Brownian particle cannot reach the 
normal diffusion regime, implicit in Einstein's relation, 
before bouncing off the pore walls. As a consequence, 
the Smoluchowski approximation fails in the vicinity of 
the bottlenecks. 

Replacing 7 with 7r in the bulk quantities Dq and 
Fq yields, respectively, Dt — kT/jT and Ft — kT/A. 
These are the 7-independent rescaling factors introduced 
in Figs. [2]|3]to characterize the inertia effects of the pore 
constrictions; 

(ii) at high drives for 

7 < 7 F = y/mF/A . (9) 
As pointed out in Ref. [T7], the large-drive regime sets 
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FIG. 3: (Color online) Rescaled diffusivity, D/Dq, versus 
F/Ft (main panel) and D/Dt versus F/Fo (inset) in the 
corrugated channel of Eq. Q with r = 1, A/yL = 0.1, 
and different 7. The scaling parameters introduced here are 
Dt = kT/jr and Fo = r y\/kT/m. The solid line in the inset 
is the heuristic power law of Eq. ( |17[ ). 



on when the length scale of the longitudinal particle dis- 
tribution grows smaller than the pore size, namely for 
F ^> Ft- In the presence a strong dc drive, the condi- 
tion 7 3> 7t does not suffices to ensure normal diffusion: 
the additional condition that A 3> If is required. Here, 
lp = mF/j 2 represents the ballistic length of a driven- 
damped particle, that is an estimate of the bouncing am- 
plitude of a driven particle against the bottleneck. Upon 
increasing F at constant 7, lp eventually grows larger 
than A and inertia comes into play. This mechanism is 
clearly responsible for the abruptly increasing branches 
of D{F) in Fig. [3} A synoptic comparison of all charac- 
teristic scaling parameters of the system is displayed in 
TableD 

In conclusion, the low- and large-drive limits are quan- 
titatively defined as f < Ft and F Ft, respectively. 
As jf was introduced to characterize the large drive (bal- 
listic) regime, clearly jp > It- This means that applying 
a large external drive makes the effects of inertia stronger. 
On the other hand, if we decrease A, while keeping F 
constant, inertia effects are controlled by 7^ rather than 
by jp, as eventually 77- > 7f- The smooth crossover 
between these two regimes is responsible for the weak A 
dependence of the fitting exponent a in Fig. [2j 

An analytical derivation of the transport quantifiers 
in the presence of strong inertial effects (low 7 and/or 
large F) proved a difficult task. This is the case, for 
instance, of the universal mobility curve in the inset of 
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FIG. 4: (Color online) Transport in a compartmentalized 
channel with r = xl/vl = 1: rescaled mobility, 7/1, versus 
F/Fo (a) and versus 7/7T (b); diffusivity, D/Dq, versus F/Ft 
(c). The remaining simulation parameters are reported in the 
legends. The relevant scaling parameters are Fo = jy^kT/m, 
Ft = kT/A, and 7t = VmkT / A. Inset of (a): jfi versus 
F/Ft for different r and A/yL- Inset of (b): 7/i versus 7/7T 
for different A/yL and r. Inset of (c): D/Do versus A/yL for 
F/Ft = 2 10 3 . The dotted curves represent the approximate 
analytical expressions of Eqs. (Ill and (15 1, respectively for 



the mobility in (b) and the diffusivity in (c) (main panel and 
insets). In (b) the quantity 7/i|oo was estimated from the hor- 
izontal asymptotes. Note that 7/^00 (see Sec. |TlJ is known to 
be proportional to A for F — > 00 [horizontal arrows, Eq. Q], 
and to llnAI for fm. 
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Fig. [2ja) . To gain a deeper insight on this and related 
issues we address next the particular case of a rectangular 
compartmentalized channel. 



IV. SEPTATE CHANNELS 

The role of inertia in compartmentalized channels is il- 
lustrated by the plots of Fig. |4j In panel (a) the rescaled 
mobility curve "/11(F) at low damping exhibits a horizon- 
tal asymptote for F — > 00. However, in comparison with 
the overdamped case reported in Sec. [TTJ such an asymp- 
tote is proportional to A only for relatively narrow pores 
(see also inset) and is strongly suppressed with decreasing 
7. The dependence of the mobility on the damping con- 
stant is better illustrated in panel (b), where 7// linearly 
increases with 7 before reaching the limit predicted in the 
overdamped regime [14H17) . Similar behaviors were ob- 
served both at low (inset) and large drives (main panel). 
For large drives the rescaled mobility actually converges 
toward the estimate 7A t |oo 

of Eq. @. 

The dependence of the rescaled mobility on the com- 
partment geometry is further illustrated in the inset of 
Fig. Qb): in the zero-drive limit and for low 7, the mo- 
bility is proportional to the aspect ratio r = xl/ul and 
the pore size, A, as already reported for the corrugated 
channels of Fig[2j[d). 

Contrary to the smoothly corrugated channels of 



Sec. Ill the drive dependence of the diffusivity is appar- 
ently not much affected by inertia. As shown in Fig.|4]jc), 
the curves D(F) keep diverging quadratically with F, 
irrespective of the compartment size and the damping 
constant, like in the overdamped limit. In the notation 
introduced above for the corrugated channel, D(F) scales 
like F@ but, contrary to Fig. [3J here /? = 2. The power- 
law dependence of D(F) on the pore size and the channel 
width is displayed in the insets of Fig. |4|c). 

Although of lesser applicability, septate channels have 
a practical advantage over smoothly corrugated channels, 
as they are characterized by distinct time scales, which 
often allow convenient analytical approximations. The 
problem under study is no exception. 

Let us consider first the rescaled mobility at large 
As anticipated in Sees. [T] and III 



us 

drives, F > F T 



two 



time scales control the particle current through the chan- 
nel: (i) the bulk relaxation time, m/7, and (ii) the bal- 
listic time across the pore, m/jx- The latter is a mea- 
sure of transient effects that may be detected only at the 
shortest distances, here, the pore width. To bridge the 
above time scales we introduce the effective relaxation 
time r e ff = rn/7 c ff, where the effective damping constant 
is defined as 



7cff = 7(1 + 7t/t)- 



(10) 



Correspondingly, the rescaled mobility function can be 
approximated to 



7/i 



7/4 



where 7/^00 denotes the rescaled mobility in the over- 
damped limit, 7 —¥ 00. Despite its being a simple inter- 
polating formula, Eq. ( 11 ) fits quite closely the simulation 
curves of Figs. 4^a) and (b) for large drives (main panel). 
Note that the horizontal asymptotes for large 7 coincide 
with the expected values of 7/i|oo> whose dependence on 
the compartment geometry, noise and drive intensity is 
analytically known |15H18j . 

Let us consider next the rescaled mobility at low drives, 
F <C Ft- For F = the transport quantifiers 7/Z0 and 
D(0) can be formally expressed in terms of the mean exit 
time, f e , of the Brownian particle out of a single compart- 
ment, namely, D(0) — x\/Af e and /io = D(Q)/kT, see 
Eq. Q. An analytical expression for f e as a function of 
the compartment geometry is only available in the over- 
damped dynamics approximation [27] . In the absence 
of a fully analytical treatment, we interpret the numer- 
ical results shown in the inset of Fig. |4jb) by assuming 
a ID collisional dynamics along the x-axis. At very low 
damping and F = 0, the particle bounces off the same 
compartment wall with rate 2v t h/xL (attack frequency), 
but only a fraction A/y^ of such collisions leads to a pore 
crossing. As a consequence, f e ~ a;Lyi/2Ay / fcT/m and 



7Mo 



1X L 



A 



r 7 

2VmkT Vl 2 7 T ' 



(12) 



which reproduces with the linear fit in the inset of 
Fig. |4^b) . Note that such a qualitative argument applies 
to the weakly corrugated channels of Fig. 2|c), as well. 
In that case, however, vth must be replaced by (2/7r)v t h, 
to account for an almost isotropic 2D distribution of the 
ballistic trajectories inside the compartment; correspond- 
ingly, the factor 1/2 on the r.h.s. of Eq. (12) should be 
changed to tt/4, see Fig. [2jd). 

The scaling law of the diffusivity at large drives, 
D(F) oc F@ with P = 2, can be quantitatively deter- 
mined by generalizing an argument originally introduced 
for the overdamped regime |17j . For large F, the instan- 
taneous particle velocity, v(t) = x(t), switches between a 
locked mode with vq = 0, as it sticks against a compart- 
ment wall, and a running mode with = F/j, as it runs 
along the central lane of the channel. In view of Eq. Q , 
it is clear that the particle spends a fraction 1 — A/it/j, 
of the time in the locked mode, and the remaining A/t/j, 
of the time in the running mode. The random variable 
v(t) can thus be modeled as a dichotomic process with 
subtracted autocorrelation function HE\ 



C(t) ee lim[(v(t + s))-(v)}[(v(s))-(v)] 

ft 

TqToo _ 



(woo - v ) exp 



(ii) 



where f = (1 — A/y L )f and = (A/y L )f are the av- 
erage permanence times, respectively, in the locked and 
running mode; their sum, f, is the relaxation time con- 
stant of the dichotomic process, still to be determined. 
The spatial diffusivity D(F) can be obtained by integrat- 
ing C(t) over time and then making use of the explicit 
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expressions for vq, i>oo, to, and Too, namely 



D(F) = / C(t)dt = 
Jo 



A/ A 

Vl V 2/i 



(13) 



To determine the unknown time constant f = (jji,/ A)f oc , 
we notice that a particle remains in the running mode for 
a time of the order of the time it takes to diffuse out 
of the central channel lane, namely, for low damping, 

2D Q r O0 ^- 4 [{y L + Af-(y L -/\f]. (14) 

By inserting the analytical expression for Too thus derived 
into Eq. ( 13 ) and taking for simplicity the limit of narrow 



pores, A<?/i, one arrives at 



D{F) 



A 



F 

JPrp 



(15) 



This expression is independent of 7 and well reproduces 
all simulation data of Fig. |4|b) at large F or, more pre- 
cisely, under the condition that 7 -C Jf ■ 

On comparing the asymptotic laws for the diffusivity 
at 7 — y 0, Eq. (151, and at 7 — > 00, Eq.(|6|, one would 
expect D{"f — > 0)/D(j —> 00) = A/j/l- On the contrary, 
in Fig. [4^c) we immediately notice that all D(F) curves 
approach the same asymptotic scaling law, Eq.([6]). As 
discussed for the corrugated channels, the overdamped 
diffusion scaling law, Eq. (15), holds only under the con- 
dition that 7 > 7j?. Correspondingly, here as well, in- 
creasing F such that F > 7 2 A/m, or lp ^ A, makes 
inertia effects on confined diffusion emerge (though in a 
less dramatic way). 



V. DISCUSSION 

The comparison between transport properties in 
smoothly and sharply corrugated channels is suggestive 
of the role played by the channels profile in the pres- 
ence of inertia. In principle, both channel geometries of 



Sec. Ill and IV could be reproduced by means of one 
parametric profile function, say, 



w v (x) 



A 



(16) 



with tunable exponent 77 |29j . This function coincides 
with w(x), Eq. ([5]), for 77 = 2 and approaches a rect- 
angular compartment ii x for 77 — >• 0. The divide 
between smoothly and sharply corrugated channels can 
be set at 77 = 1, where the two sides of the bottleneck 
profile change from concave, for 77 > 1, to convex, for 
77 < 1. Such change in the pore geometry affects, for in- 
stance, the drive dependence of the rescaled mobility at 
low damping, see Fig. [5} All curves 7/i(-F) with 77 > 1 
decay with the same approximate power law as reported 
in Fig. ^b) for 7/ = 2. For 77 < 1, instead, the behavior 
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FIG. 5: (Color online) Channel with tunable corrugation, 
Eq. (16 1: rescaled mobility, 7^1 versus F/Ft for different 
77. Other simulation parameters: r = 1, A/yL = 0.1, and 
7/7T = 0.8. The dashed line is the power-law (F/F T )~ a/2 , 
with a — 1.4 drawn in Fig. |2|b) for rj — 2. 



of 7/j(F 1 ) is as in Fig. [2^b) (sinusoidally corrugated chan- 
nel), at low F, and in Fig.Qa) (septate channel), at large 
F. Without further analyzing the 77 dependence of the 
transport quantifiers, we now discuss certain differences 
and similarities between sinusoidal and septate channels. 

In both types of channels, the diffusivity grows asymp- 
totically with the drive according to a power law, D{F) cx 
F 13 , where j3 — 2 for septate channels and j3 = 1 for sinu- 
soidal channels. For smoothly corrugated channels this 
result may come as a surprise, since, for large damp- 
ing, the diffusivity is known to approach its bulk value, 
D(F —> 00) = D(0). The asymptotic power law with 
/3 = 1 illustrated in the inset of Fig. [3] results from the 
enhanced trapping effect due to the underdamped par- 
ticle bouncing back and forth against the compartment 
walls. A simple qualitative argument yields |26j 



D(F) ^ nx L F ^^Vl F 
D T ~ 8 kT r 8 A F T 



(17) 



in rather good agreement with our simulation data (see 
inset of Fig. [3]). Note that, lowering the temperature, for 
small damping D(F) diverges like T" 1 / 2 , which means 
that diffusion is the result of chaotic ballistic collisions. 

Inertia corrections to the drive dependence of the diffu- 
sivity are not as dramatic in septate channels, Fig. |4^c) , 
as they appear in corrugated channels, Fig. [3j This ex- 
plains why the role of the threshold 7f, Eq. ([9]), is less 
prominent for sharply corrugated channels. On a closer 
look, however, one notices that, on increasing 7, the 



data points for D(F)/Dp approach the predicted oblique 
asymptote of Eq. ( 15 1 at larger and larger F, consistently 



with the large-drive inertial regime 7 < 7j?. 

Finally, we notice that septate and sinusoidally corru- 
gated channels also differ in the large-drive behavior of 
their mobilities at low damping. While in septate chan- 
nels 7/x(i ? — ¥ 00) was shown to approach a small but 
finite value, 7/^|oo(7/7t), see Eq. (If), the mobility in 



sinusoidal channels was numerically fitted by the scal- 
ing law 7/i(-F — > 00) ~ (7 '/7f) q : , where a is an increasing 
function of A with a(A -> 0) = 2 [Fig.^c)]. This means 
that in sinusoidal channels (v(F)) [and not 7/j,(F)] tends 
to a finite asymptotic value. By definition, the rescaled 
mobility can be written as 7/i ~ (xr,/fd)('j/F), where fd 
denotes the mean drift time of a particle across a com- 
partment in the presence of a strong drive. Accordingly, 
as a — ¥ 2, the drift time fd becomes insensitive to the 
(large) drive, which hints at an emerging ballistic dy- 
namics [3U]. We also remark that the above scaling law 
for 7/! (-F -> 00) applies to all pore geometries with 77 > 1 
[29] . see Fig. [5j for septate channels, 77 — > 0, such a scal- 
ing law, with f < a < 2, closely reproduces the decaying 
branch of the mobility curves displayed in Fig. |4|a). 



VI. CONCLUSIONS 

The main result of this work is that for real physical 
suspensions flowing through confined geometries, both in 
biological and artificial systems, pore crossings grow in- 
creasingly sensitive to the suspension fluid viscosity with 
decreasing the pore radius. With respect to previous at- 
tempts at incorporating finite-mass effects in the analysis 



of Brownian transport through corrugated narrow chan- 
nels [3TH33] . we stress that the inertial effects reported 
here are not of mere academic interest |26| . 

Inertial effects can be directly observed, for instance, 
in a dilute solution of colloidal particles driven across a 
porous membrane or an artificial sieve [Ml ISS] . On the 
other hand, channel profiles at the micro- and nanoscales 
can be tailored as most convenient |36j . As detailed in 
Ref. the experimental demonstration of inertial ef- 
fects on Brownian transport through narrow pores is to 
become accessible when manipulating artificial particles 
of micrometric size by means of well established exper- 
imental techniques [3"T1 - |4"0"] . For nano-particles, like bio- 
logical molecules, detecting such effects will require more 
refined experimental setups. 
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